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ABSTRACT 

In this paper, we use the general theory worked out within the past few 
years for the structure and the evolution of low-mass stars to derive the stellar 
mass-function in the Galactic disk down to the vicinity of the hydrogen-burning 
limit, from the observed nearby luminosity functions. The accuracy of the 
mass-magnitude relationships derived from the afore-mentioned theory is 
examined by comparison with recent, accurate observational relationships in 
the M-dwarf domain. The mass function is shown to flatten out below ~ 1 M Q 
but to keep rising down to the bottom of the main sequence. Combining 
the present determination below 1 M Q and Scalo's (1986) mass function for 
larger masses, we show that the mass function is well described over the entire 
stellar mass range, from ~ 100 M Q to ~ 0.1 M Q , by three functional forms, 
namely a two-segment power-law, a log-normal form or an exponential form, all 
normalized to the Hipparcos sample at 0.8 M & . 

Integration of this mass function yields a reasonably accurate census of 
the entire stellar population in the Galactic disk, and its volume and surface 
mass-density 

Subject headings: Dark matter, low-mass stars, luminosity function, mass 
function 
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1. Introduction 

A correct determination of the Galactic mass function (MF) over the entire stellar 
domain is essential for several reasons. First it yields the complete census of the stellar 
population in the Galaxy and thus its contribution to the Galactic baryonic content and 
mass-to-light ratio. Second, the normalization and the slope of the stellar MF near the 
hydrogen-burning limit provides the boundary condition to infer the brown dwarf MF, 
assuming the MF is continuous at the H-burning limit. Third, the determination of 
the global stellar plus substellar MF provides an essential diagnostic to understand the 
formation of star-like objects. 

Star counts, which a few years ago still relied essentially on giants and sun-like stars, 
i.e. objects with mass m>lM & , now include information on the M-dwarf population 
in the Galaxy. As discussed below, modern observations now probe the M-dwarf stellar 
distribution down to the bottom of the main sequence (MS), close to the hydrogen-burning 
limit. Moreover, we know that this distribution extends into the substellar regime since 
over a hundred field brown dwarfs have now been discovered in the solar neighborhood. 

In the present paper, we define the initial mass function (IMF) £ as the number of 
stars N, or the stellar number- density per cubic parsec n, formed initially per mass interval 
[to, m + dm] : 



e( \ dN 
^ = d^ 



The MF can also be described as the number or the density of stars per interval of 
logm, as defined originally by Salpeter (1955), with the straightforward relation: 
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flog(m) = - = (LnlO) mf(m) (2) 

aiogm 

The quest for the determination of the stellar IMF takes root in the seminal paper by 
Salpeter (1955) and in the following major contributions by Miller & Scalo (1979), Scalo 
(1986), and Kroupa, Tout & Gilmore (1991, 1993). All these papers agree for an IMF 
reasonably well described by a power-law form £(m) oc m~ a with an exponent a ~ 2.35 for 
stars with m> 1 M , the so-called Salpeter IMF, although the more detailed analysis of 
Scalo (1986) suggests a steeper slope, with a ~ 2.7. A recent analysis by Kroupa (2001) 
suggests a ~ 2.3 but, as noted by this author, the previous Scalo value is recovered if 
one assumes a standard fraction of observationally unresolved binaries for these massive 
stars. On the other hand, the IMF seems to flatten significantly at the low-mass end, as 
suggested initially by Miller & Scalo (1979). Indeed, the detailed analysis of Kroupa et al. 
(1993) suggests a significant flattening of the IMF below ~ 0.5 M Q , with a ~ 1.3, although 
the question of whether the IMF keeps rising down to the H-burning limit or decreases 
below a certain mass around ~ 0.2 M & remains unclear. In fact the very functional form 
of the IMF in the M-dwarf range, power-law, log-normal or other form, is presently still 
undetermined. It is important to note, in the present context, that the approach of Kroupa 
and collaborators, who first addressed this problem in the M-dwarf range, is essentially 
empirical, in the sense that these authors used an empirical mass-luminosity relationship 
derived from observations of nearby binaries (Popper, 1980). 

The aim of the present paper is to demonstrate those results, in deriving a stellar 
MF from a consistent stellar evolution theory, and to explore the possibility of various 
functional forms for the IMF in the low-mass range. As mentioned below, this theory 
describes accurately the very mechanical and thermal properties of M-dwarf stars, and 
yields reliable mass-magnitude relationships, allowing age and metallicity determinations. 
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The present paper presents the derivation of the stellar IMF for the Galactic disk over the 
entire stellar mass range, with particular emphasis on the low-mass part of the distribution 
(m < 1 Mq). The accuracy of the stellar models, more particularly of the mass- magnitude 
relationships, is examined in §2 by comparison with the most recent observationally 
determined relationships in different optical and infrared passbands. In §3, we examine the 
different M-dwarf luminosity functions presently available for probing the disc low-mass 
star distribution. The mass functions are derived in §4. In this section, we carefully 
examine the uncertainties due to the mass-magnitude relationship and we derive analytical 
parametrizations of the MF. The inferred stellar mass budget is determined in §5 while 
section 6 is devoted to the conclusion. The extension of these calculations into the brown 
dwarf regime, where age becomes an extra degree of freedom besides mass, will be presented 
in a subsequent paper. 

2. The mass-magnitude relationships 

The key ingredient in the determination of the MF is the mass-magnitude relationship 
(MMR). As too rarely stressed, the MMR is the cornerstone to transform the observable 
quantity, the luminosity function (LF) = dN/dM, i.e. the number of stars N per absolute 
magnitude interval [M,M + dM], into a MF. The derivation of observable MMRs for 
M-dwarfs is a formidable task for several reasons. First of all, we need to determine the 
absolute magnitudes of the objects, which implies accurate parallaxes, and thus nearby 
distances. Second, M-dwarfs are by definition intrinsically faint, in particular near the 
bottom of the MS, where the determination of the MMR is most crucial since the luminosity 
drops by orders of magnitude with the mass when approaching the hydrogen-burning 
limit. At last, the determination of the mass requires binary systems, which reduces the 
already weak statistics by about a factor 2, and dynamical information, i.e. a long enough 
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time basis. All these difficulties render the determination of an observed M-dwarf MMR 
a formidable challenge for astronomers. A first compilation of mass-luminosity data in 
the M-dwarf domain was published by Popper (1980), and was subsequently extended by 
Henry & McCarthy (1993). These authors used speckle interferometry to obtain MMRs 
in the V-, J-, H- and K-bands. The determination of the V-magnitude of these objects 
has been improved recently by using the HST (Henry et al., 1999), reducing appreciably 
the uncertainty in the m-My relation. The observationally-determined MMRs were fitted 
by simple polynomial expressions. As mentioned by the authors themselves, these fits 
are to be taken with great caution and must be used as guidelines only, for they do not 
take into account age and metallicity-dispersion, characteristic of the young disk/old disk 
population. The Henry & McCarthy (1993) sample has been improved significantly recently 
by Delfosse et al. (2000). These authors combined adaptative optics images and accurate 
radial velocities to determine the mass of about 20 objects between ~ 0.6 and ~ 0.09 M 
in the V-, J-, H- and if-bands. The masses of the Delfosse et al. sample are obtained with 
accuracies of 0.2 to 5% and the distance is determined also with high precision. 

Although such data provides extremely useful empirical MMRs, the MMR is generally 
derived from theoretical models aimed at describing as accurately as possible the structure 
and the evolution of the objects under study. Indeed the shape of the MMR reflects the 
very mechanical and thermal properties of these objects (see Chabrier & Baraffe, 2000 for a 
recent review). The accuracy of the stellar models and of the derived theoretical MMRs is 
thus a crucial issue since the transformation of the LF into the MF involves the derivative 
of the MMR: 



dN dN dm 1 

dm~ {m) = { dM^m~y X [ m^m) ] (3) 

where N is the number of stars, m the mass and M\ the absolute magnitude in a given 
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passband. 

For a claimed MF determination to retain any degree of reliability, it must rely on 
MMRs, and thus on stellar models, which have been demonstrated to accurately reproduce 
observational data in all available observational diagrams, color-color, color-magnitude, 
mass-spectral type and most importantly mass-magnitude. This is mandatory for any 
MF determination. Only recently has a generation of low-mass stars (LMS) and brown 
dwarfs (BD) evolutionary models emerged which relies on consistent calculations of the 
internal and atmospheric properties of these objects, with no adjustable parameter^ and 
thus provide consistent mass-age-magnitude-color relationships. The physics of these 
models and extensive comparisons with observations in different observational diagrams 
have been discussed in various papers (Chabrier & Baraffe 1997; Baraffe et al 1997, 
1998; Chabrier et al., 2000; Chabrier & Baraffe, 2000 and references therein). Although 
improvement is still needed for a complete agreement with all the observed magnitude-color 
diagrams (see Baraffe et al., 1998, hereafter BCAH98 and discussion in §4.1), these 
models improve significantly the situation, in particular in the infrared, and reach now 
quantitative agreement with the observations. Figure 3 of Delfosse et al. (2000) shows a 
comparison of the MMRs issued from the afore-mentioned general theory with the Delfosse 
et al. (2000) recent data in the V, J, H and i^-bands. The agreement in the J, H and 
i^-bands is essentially perfect (< la). The situation is not as satisfactory in the V-band, 
with a systematic offset between theory and observation below ~ 0.3 M , My > 12. The 
consequence of such an uncertainty on the derivation of the MF will be examined carefully 
in 84.1. 



1 Recall that below ~ 0.7 M the objects are dominantly or even entirely convective so that 
variations of the mixing length parameter is inconsequential on the evolution (see Chabrier 
& Baraffe, 1997 for a detailed discussion) 



This general excellent agreement between theory and observation brings confidence in 
the description of the mechanical and thermal properties of low-mass stars down to the 
vicinity of the H-burning limit. It allows a more detailed study of the effects of age and 
metallicity on the MMR. 

• age effect : objects with m>0.1 M reach the MS after ~ 0.7 Gyr (see e.g. Table 1 
of Chabrier & Baraffe, 2000) and, for m<O.8M , stay there for a Hubble time, so that 
their position is fixed in the mass-luminosity diagram. Younger objects in this mass-range 
will still be on their pre-MS contraction phase and will be more luminous than main 
sequence objects with the same mass (see Figure 3 of BCAH98). Assuming a constant 
stellar formation rate and an age for the Galactic disk to ~ 10 Gyr, the fraction of such 
young objects in the local population, however, should represent at most a few percents. 

• metallicity effect : The Hipparcos color-magnitude diagram indicates that ~ 90% of 
disk stars have abundances within ±0.2 dex of the solar value (Reid, 1999) so the spread 
of metallicity in the solar neighborhood should not affect significantly the derivation of 
the MF through the MMR. The consequence of variations of the metallicity over a range 
A[M/H] = ±0.5 on the MF is examined in Mera, Chabrier & Baraffe (1996). A lower 
metallicity yields a steeper slope since the lower the metallicity the brighter the star for a 
given mass or, conversely, the lower the mass for a given magnitude. As mentioned by these 
authors, however, variations of the MF obtained within the afore-mentioned metallicity 
range remain within the Poisson error bars of the value determined for a solar abundance. 

As discussed in detail in BCAH98 (§3) and Delfosse et al. (2000), the MMR in 
infrared bands is very weakly affected by metallicity variations in the afore-mentioned 
range. Indeed, metallicity effects for a given mass cancel out with the effect of decreasing 
effective temperature on the spectral energy distribution in infrared bands whereas they 
add up in the V-band. Therefore, the derivation of the MF from a LF determined directly 
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in near-IR bands should be insensitive to biases due to a metallicity spread in the sample. 
Unfortunately there is presently no direct determination of the LF in such bands. The 
K-band LF described in §3 is derived from the V-band LF through various color-color 
relations and is thus hampered by uncertainties du to these transformations. The thick-disk 
population ([M/H] ~ —0.5), with a scale height h ~ 0.7 kpc and a local normalization of 
~ 5% (Haywood, Robin & Creze, 1997) should make a negligible contribution to the local 
sample. 

• grain formation : below T e ff m 2800 K, grain condensation in the atmosphere of 
low-mass objects affects drastically their spectral distribution, and thus the mass-magnitude 
relationship (Chabrier et al., 2000). This corresponds to objects less massive than 0.1 M & 
and is not consequential for the present study. 

3. The Galactic disc M-dwarf luminosity- function 

The determination of the M-dwarf LF in the Galactic disk down to the bottom of 
the MS has been the source of particular activity within the past recent years. The 
determination of the absolute LF requires the determination of the distance of these objects. 
The easiest way to determine the distance is by knowing the trigonometric parallax of 
the object, which implies a search within near distances from the Sun, typically d < 20 
pc for the bright part of the LF (My < 9.5), a few parsecs for the faint end. This yields 
the so-called nearby LF. The main caveat of the nearby LF is that, given the limited 
distance, it covers only a limited volume and thus a limited sample of objects. This yields 
important statistical indeterminations at large magnitudes [My > 12). On the other hand a 
fundamental advantage of the nearby LF, besides the reduced error on the distance, i.e. on 
the magnitude, is the accurate identification of binary systems. Other determinations of the 
disk LF are based on photographic surveys, which extend to d ~ 100 — 200 pc from the Sun, 
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and thus encompass a significantly larger amount of stars. However, photometric LFs suffer 
in general from significant Malmquist bias and, most importantly, the low spatial resolution 
of photographic surveys does not allow the resolution of binaries at faint magnitudes. An 
extensive analysis of the different nearby and photometric LFs has been conducted by 
Kroupa (Kroupa, 1995; see also Reid & Gizis, 1997, and references therein). As shown by 
this author, most of the discrepancy between photometric and nearby LFs for My > 12 
results from Malmquist bias and unresolved binary systems in the low-spatial resolution 
photographic surveys. 

Kroupa (2001) derived a nearby LF, $ near , by combining Hipparcos parallax data, 
which is essentialy complete for My < 12 at r=10 pc, and the sample of nearby stars 
(Dahn et al., 1986) with ground-based parallaxes for My > 12 to a completeness distance 
r=5.2 pc. Reid and Gizis (1997, RG) extended this determination to a larger volume and 
determined a nearby LF based on a volume sample within 8 pc. Their local sample is 
drawn from the most recent version of the Catalogue of Nearby Stars (Gliese & Jahreiss, 
1991). Most of the stars in this survey have parallaxes. For late K and M dwarfs, however, 
trigonometric parallaxes are not always available; in such cases, these authors use a 
spectroscopic TiO-index vs My t0 relation to estimate the distance (Reid et al., 1995). This 
sample was revised recently with Hipparcos measurements and new binary detections in the 
solar neighborhood (Delfosse et al., 1999) and leads to a revised northern 8-pc catalogue 
and nearby LF. RG argue that their sample should be essentially complete for My < 14, 
with all stellar companions resolved. About 35% of the systems in the RG sample are 
multiple and ~ 45% of all stars have a companion in binary or multiple systems. 

The HST LF (Gould, Bahcall & Flynn, 1997), §hst, extends previous photometric 
surveys to an apparent magnitude I < 24. The Malmquist bias is negligible because all stars 
down to ~ 0.1 M & are seen through to the edge of the thick disk. A major caveat of any 
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photometric LF, however, is that the determination of the distance relies on a photometric 
determination from a color-magnitude diagram. In principle this requires the determination 
of the metallicity of the stars since colors depend on metallicity (see e.g. Figure 12 of 
Chabrier & Baraffe, 2000). This is particularly true for the HST LF since a substantial 
number of stars in the sample belong to the thick disk and have metal-depleted abundances. 
Moreover, as mentioned above, the HST small field of view does not allow binary resolution 
so that the HST misses essentially all the binaries, yielding only the determination of the 
stellar system LF. Although, as mentioned above, unresolved binary systems help resolving 
the differences between the nearby LF and the HST LF, they cannot bring the two LFs 
into agreement and at least a factor 2 difference remains at the faint end {My > 14) of the 
LF and of the MF (see Mera, Chabrier & Schaeffer, 1998, their figure 1, with Figures [1] and 
H below). The source of this disagreement is presently unclear and will be the object of a 
further study. 

The illustration of the disagreement between the three afore-mentioned LFs, $5.2 pc , &8pc 
and &hst in the V-band at faint magnitude {My > 12) is shown in Figure 15 of Chabrier & 
Baraffe (2000). Whereas the photometric LF exhibits a pronounced peak around My « 12 
[Mj p« 10), with a sharp drop below, the two nearby LFs, although still peaking at this 
magnitude, remain rather flatfish up to My ~ 15, predicting a substantially larger number 
of stars in this region, a disagreement which cannot be explained only by unresolved 
binaries, as mentioned above. The reasonable agreement between the two nearby LFs, $5.2 pc 
and $8pc down to My ~ 15, about the limit of completeness claimed for the 8-pc sample, 
brings some confidence in these determinations. The disagreement between these two LFs 
at fainter magnitude stems very likely from incompleteness of the 8-pc sample. As shown by 
Henry et al. (1997; their figure 1) and Delfosse et al. (1999), the M-dwarf sample becomes 
substantially incomplete for distances larger than 5 pc. Second, several M-dwarfs in the RG 
sample, which represent the faint end population of the LF, lack parallax determination 
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and have distances estimated from the afore-mentioned spectroscopic relation. This relation 
can be uncertain by more than 1 magnitude for My > 10 (see Figure 3 of Reid et al., 
1995), a major source of Malmquist bias. On the other hand, the last bins in the $5.2 pc 
LF {My > 16), which include only 5 stars or less each, might be contaminated by one or 
several young brown dwarfs, by young pre-MS objects, or by low-metallicity (brighter) 
very-low-mass stars. It is unfortunately not possible at the present stage to resolve this 
issue. We have thus calculated MFs from both nearby LFs, which represent the most 
complete LFs for M-dwarfs presently available. 

RG have also derived the LF in the K-band. However, as mentioned above, there is 
no direct determination of the K-band LF, at least for M-dwarfs, and these authors used 
(V — I) vs (V — K) or BCj vs (V — I) relations, or M^ vs My relations for stars with only 
an absolute magnitude, fitted to the Leggett (1992) data (see below). 

4. The mass function 

Figure [l] displays the MF £(logm) calculated with the theoretical MMRs mentioned 
in §2 from the V-band LFs $y 52pc (solid line), ^v Spc (dash line) and from the V-band LF 
converted in K $^ 8 c (dotted line), in 10 3 pc 3 . The error bars are the ones quoted by the 
authors for the RG revised sample, and are la Poisson errors for the 5.2 pc sample. Note 
that the MF from the RG 8-pc sample has been calculated from a star-by-star analysis, 
from the sample kindly provided by Neill Reid, not from eqn.(l). We verified, however, that 
both methods yield similar results, within the error bars. The first interesting result is the 
fact that the MFs derived from the two different V-band LFs agree very well, within la in 
average, down to logm = —0.8 (m = 0.15 M Q ). The highest mass bin of the 8-pc MF lies 
several a's above the MF derived from the Hipparcos sample and must be considered with 
caution. The 2-a discrepancy between these two MFs below logm ~ —0.9 (m = 0.12 M , 
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My ~ 14.5 at 5 Gyr) stems very likely from the acknowledged incompleteness of $g pc at 
these magnitudes. On the other hand, the last bin of the 5.2-pc MF may be contaminated 
by a massive BD or a young pre-MS star, as mentioned previously, and thus should be 
considered with limited confidence. The MF derived from the K-band $8 pc LF differs at the 
2<t level from the two previous ones, and exhibits a more flatfish general behaviour. There 
is a priori no reason for such a difference since the sample is exactly the same as for the 
$y 8 c LF. There are two possible explanations for such a behaviour: 

i) the disagreement might stem from discrepancies in the m-My behaviour, as 
mentioned above. This will be examined in detail in the next subsection. 

ii) only 68 of the 106 systems in the 8-pc sample have wide-field K-band imaging. For 
objects in the sample lacking IR observations, the Mk absolute magnitude is determined 
either from a (V — K) versus (V — I) relation or from a Mk versus My relations, for those 
lacking V,R,I photometry, based on the Leggett (1992) red dwarf sample (see RG for 
details). Although such relations provide reasonable average transformations, they do not 
take into account the color spread in the observational diagram. In particular the ~ 0.2-0.3 
mag spread in color in the observed (V — K)/(V — I) diagram for 2.5 <V — J<3.5, i.e. 
0.4 >m/M & > 0.2. This translates into a significant uncertainty in the derived K-band LF. 
In order to examine this effect, we have recalculated the K- magnitude in the RG sample for 
objects with V, R, I photometry with a slightly different relation: 



(V -K) = l.075 + lAx(V-I) (4) 

which seems to provide a better mean fit of the observed (V — K)/(V — I) relation than the 
one used by RG. The MF derived from this new $/^ g c is shown in the next subsection. 
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4.1. Uncertainty due to the mass-magnitude relation 

As stressed previously, the reliability of the MF depends crucially on the accuracy 
of the MMR. As mentioned in §2 and discussed in BCAH98 and Delfosse et al. (2000), 
there is a discrepancy between the theoretical mass-My relationship derived by BCAH98 
and the observationally derived relation below ~ 0.35 M , the models being systematically 
brighter than the observations, with a maximum ~ 0.5 magnitude discrepancy around 
m = 0.2 Mq. This was anticipated from an identified shortcoming of the present theory 
in the V-band (see BCAH98), mainly due to still incomplete molecular linelists and thus 
opacity coefficients in the optical. Therefore, the models underestimate the mass for a 
given magnitude below ~ 0.35 M . The effect, however, remains modest: for example a 
magnitude My = 12 corresponds to an observed mass m cz 0.28 M Q , whereas the theory 
predicts m ~ 0.25 M ; the maximum disagreement around My = 13 corresponds to an 
observed mass m ~ 0.2 M , whereas the theory predicts m ~ 0.17 Mq. In order to quantify 
the effect on the MF, we have recalculated it with the parametrized mass- My relation of 
Delfosse et al. (2000) which provides a good mean fit of their data. This polynomial fit 
reads (see Delfosse et al., 2000): 



log(m/M ) = 10~ 3 x [0.3 + 1.87 x My + 7.614 x M v - 1.698 x M v + 0.060958 x M v ] (5) 

for M v G [9, 17] 

Figure ^| displays the MF derived from $y 5 2pc and $y gpc when using this relation, to 
be compared with Figure 0. As anticipated from the discussion above, the main effect 
of the Delfosse et al. (2000) m — My relation is to increase the number of stars in the 
mass range m ~ 0.2 — 0.4 M by slightly depressing both the higher mass tail and the 
very-low-mass tail of the distribution, yielding a slightly more flatfish distribution below 
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logm fh —0.5 (m ~ 0.3 M & ). Figure |2| displays also the MF derived from the revised 
$x 8 c LF, obtained with eqn.(4), as described above. Interestingly, the three MFs seem 
now to agree reasonably well. The 1 to 2 a difference between the MF derived from the 
5.2-pc sample and the ones derived from the 8-pc sample for logm< — 0.8 stems very 
likely from (i) incompleteness of this latter at faint magnitudes, (ii) remaining uncertainties 
in the My into Mr transformation for $^ 8 c > and (iii) from the presence of BD or PMS 
objects in the last bin of $v&2 PC - Because of point (ii), and until the M-dwarf LF is 
observationally-determined in infrared filters, we will focus in the following on the MF 
derived from the two observed V-band LFs. 



4-1.1. Parametrization of the mass fonction 

Figure |3| displays the two MFs derived from the V-band LFs with three different 
functional forms superposed, namely: 



- a power-law form (solid line; IMF1): 



nil 

&) = ^ = ^« (6) 

dm 



- a log-normal form (dot-dash line; IMF2): 



c n \ dn a i (logm - log mo) 2 . , . 

£ log m) = — = Aexp{ — } (7) 

a log m 2a z 



an exponential form (dash line; IMF3): 



£(m) = — = Am- a exp{- — V 3 } 
dm m 
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This latter form resembles a log-normal form in the low-mass range and recovers 
asymptotically to a power-law at large m, £,{jn) m>>mQ — > m~ a . The peak value corresponds 
to m p = (-^l^rriQ. Such a functional form was proposed originally by Larson (1986) in the 
case of a bimodal star formation process in galactic disks and a remnant dominated dark 
matter component in the solar neighborhood. It has been recently advocated as the IMF 
required for a dark halo containing a substantial population of stellar remnants, accounting 
for the microlensing observations towards the LMC (Chabrier et al., 1996; Chabrier 1999). 

All three functional forms are normalized at 0.8 M & on the Hipparcos value 
(^)o. 8 M = 2.8 ± 0.2 x 1O~ 2 M J9C -3 . Stars below m ~ 0.8 M did not have time to evolve 
off the MS within ~ 10 Gyr, about the age of the Galactic disk, so the present-day MF for 
stars below 0.8 M & reflects the IMF. One of the constraints to be fulfilled by the IMF is 
the so-called continuity constraint, which requires the IMF to be continuous at the mass for 
which the main sequence lifetime ttas corresponds to the age of the disk tMS ~ ^d (Miller 
& Scalo, 1979). This mass is about ~ 0.8-0.9 M & . For larger masses, the IMF seems to be 
well described by a power-law with an exponent a ~ 2.7 (Scalo, 1986), although substantial 
uncertainties remain in the exact determination on the exponent. Therefore, the forms 
(6)- (8) must be close to such a mass distribution for masses above 1 M to fulfill the 
continuity constraint. A condition reasonably well satisfied in all three cases, as shown in 
Figure |j. Table 1 gives the values of the different parameters entering equations (6)- (8) as 
well as the corresponding average mass of star formation: 



Q^^{m)mdm 

{m) = c*7 en dm (9) 

where w, nn . is the stellar minimum mass, by definition the hydrogen-burning limit, 

~ 0.07 Mq for solar-metallicity stars (Chabrier & Baraffe, 1997), and m± su is the maximum 

mass for star formation, ~ 100 Mq. 
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As seen in Figure |3|, the three analytical forms (6)-(8) provide very reasonable fits - 
within less than 2<r - of the MF, if the last bin of the MF derived from $5.2 pc is omittedQ. 
The dotted line displays the 4-segment power- law IMF derived by Kroupa (2001). 

Figure |] displays the three MFs over the entire stellar mass-range from 60 M Q down 
to the H-burning limit, and the extension in the substellar domain. As seen in the figure, 
the exponential IMF is an excellent description of the stellar distribution over the entire 
characteristic mass range, recovering a power-law tail at large masses and the log-normal 
form below ~ 1 M . Such a functional form, proposed for the dark halo IMF, as mentioned 
above, is thus not an ad-hoc prescription but represents quite well the Galactic disk stellar 
IMF. Width and normalization are of course different for the disk and the primordial 
halo (see Chabrier, 1999), as anticipated from the very different characteristic conditions 
(virial temperature, metallicity, etc.). As seen in the figure, the power-law IMF differs 
significantly from the two other forms in the substellar regime. Brown dwarf counts should 
then enable us to distinguish between these forms. Substellar objects, however, never reach 
thermal equilibrium, which defines the MS, and the age distribution must be taken into 
account when determining the IMF from the present-day observed distribution. This issue 
will be addressed in a forthcoming paper. The log-normal and the exponential forms are 
very similar: at ~ 0.01 M & (~ 10 Jupiter mass), the predicted number of objects differs 
only by about 60%. Such a difference is unlikely to be detectable for objects in the field, 
with resolved companions. The IMF derived by Kroupa (2001) is shown for comparison by 
the dotted line. 

Figure [^ displays the same MFs below 0.6 M in a linear scale. A Salpeter IMF with 
the same normalization is shown for comparison (dotted line). Such an IMF all the way 



2 Taking this last bin into account yields a significantly steeper MF, with a ~ 1.8, as 
derived in the previous analysis of Mera, Chabrier & Baraffe (1996). 
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down to the bottom of the MS would overestimate the stellar density at 0.1 M by about a 
factor of 6. 



5. The Galactic disk stellar mass budget 

Integration of the MFs derived in §4 from 1.0 M down to the H-burning minimum 
mass 0.07 M yields the LMS (generically designed here as stars below a solar mass, 
m < 1.0 M Q ) number-density ulms and mass-density Plus'- 

rl.O /-log(l.O) 

tilms = / £(m)dm = / £(logm) dlogm 

J0. 07 Jlog(0.07) 

rl.O /-log(l.O) 

Plus — \ ^{m)mdm= / £(logm) md log m (10) 

J0. 07 Jlog(0.07) 

Table 2 gives n LMS and Plms f° r the three analytical MFs given by eqns. (6)- (8). 
Adding the more massive stars with the scale height dependent present-day MF (PDMF) 
of Scalo (1986, Table 4) yields the present-day disk main sequence stellar density. 
Adding up the single white dwarf density n^ D ~ 5.5 ± 0.8 x 10~ 3 pc~ 3 (Holberg, 
Oswalt & Sion, 2001) with an average mass (itiwd) — 0.67 M , i.e. a white dwarf 
mass density pwd — 3.7 ± 0.5 x 10 -3 M pc~ 3 , a neutron star density n^s — 10~ 3 
pc -3 (Popov et al., 2000) with a mass (ijins) — 1-4 M and a red giant contribution 
n RG ~ 0.3 x 10~ 3 pc" 3 , p RG ~ 0.6 x 10~ 3 M G pc~ 3 (Haywood et al., 1997) yields the 
present-day Galactic thin-disk total stellar number- and mass-densities n* ~ 0.13 x 10 -2 
pc -3 , p+ ~ 4.30 ± 0.3 x 10^ 2 M Q pc~ 3 . Assuming a standard double-exponential thin- 
disk+thick-disk with a common scale length / = 2.5 kpc and scale heights h = 250 
pc and h ~ 760-1000 pc, respectively, for stars with m < 1.0 M , and a Scalo (1986) 
mass-dependent scale height for more massive stars in the thin-disk, with a thick-disk local 
normalization ~ 5% (Haywood et al., 1997), we obtain the stellar contribution to the disk 
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local surface density S Qa ~ 24.0 ± 2 M Q pc~ 2 . The ~ 10% error bars on Plms correspond 
to the la Poisson uncertainty on the MF at logm = —0.9. These numbers are still affected 
by the uncertainty at the very- low- mass end of the stellar MF, due to small statistics and 
incompleteness, and by the uncertainty in the exact scale height and normalization of the 
thick disk. These uncertainties, however, should not affect the present numbers by more 
than ~ 10%. An accurate determination of the MF below m ~ 0.15 M & , i.e of the LF 
at My > 14, Mk > 9, is strongly needed to reduce this uncertainty and to nail down the 
normalization of the stellar MF at the H-burning limit. 

6. Conclusion 

We have derived the low-mass star mass-function characteristic of the disk stellar 
population from parallax-determined LFs down to the vicinity of the hydrogen-burning 
limit. This stellar MF is derived from a consistent stellar evolution theory which accurately 
describes the very mechanical and thermal properties of these stars, and yields reliable 
mass-magnitude relationships. The consequences of the remaining uncertainties in the 
theoretical mass- My relationship on the determination of the MF have been examined 
carefully. The MFs derived from both the 5.2-pc and 8-pc LFs in the V-band agree 
fairly well (< la) down to ~ 0.15 M . The MF derived from the K-band LF is in 
reasonable agreement with the previous ones, but its accuracy is hampered by photometric 
or magnitude transformations from Mk to My based on simple mean fits. The direct 
determination of the LF in the K-band would be of tremendous interest, since the dispersion 
in the mass-magnitude relation due to the metallicity spread is negligible in this band. 
Incompletude is likely to affect the 8-pc sample at magnitudes corresponding to the 
above-mentioned mass. Age effects are consequential only for objects below ~ 0.1 M & and 
younger than ~ 0.5 Gyr and thus do not affect significantly the main sequence stellar mass 
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function above this limit. This can affect, however, the very faint end of the nearby 5.2-pc 
LF, leading to the misidentification of a statistically significant number of young brown 
dwarfs or pre-MS very-low-mass stars as main sequence low-mass stars. Nevertheless, 
agreement between the various nearby surveys down to My ~ 15 brings confidence in the 
present determination of the stellar MF down to about 0.15 M & . 

The MF is found to flatten significantly below about ~ 1.0 M , compared with a 
Salpeter MF, and again below ~ 0.5 M Q , as noted originally by Miller & Scalo (1979) and 
Kroupa et al. (1993), but keeps rising slowly down to the bottom of the main sequence. 
Incompleteness and limited statistics prevent the unambiguous determination of the value 
of the mass function at the very hydrogen-burning limit. Larger statistics is needed at 
faint magnitudes, a difficult task since parallax-based surveys are likely to be affected 
by incompleteness beyond ~ 5 pc and photometric surveys lack binary resolution, not 
mentioning possible erros due to Malmquiest bias or color-magnitude transformations. 

The IMF is well described over the entire stellar mass range, i.e. four orders of 
magnitude in mass, by three different functional forms, namely a two-segment power-law, 
a log-normal and an exponential. Interestingly, this latter form, which has been advocated 
recently to suggest a substantial population of old stellar remnants in the dark halo, seems 
to offer the best compromise between the high-mass tail and the low-mass end of the stellar 
distribution. Only in the substellar regime does the power-law IMF differ significantly from 
the two other forms, predicting a substantially larger number of field brown dwarfs. This 
issue will be addressed in a forthcoming paper. The previously suggested bimodal behaviour 
for the IMF, with a maximum around ~ 0.3M Q and a drop below (Scalo, 1986), seems to 
be excluded. It is very likely an artifact due to (i) unresolved binaries in the photometric 
luminosity functions (see Kroupa et al., 1993), (ii) small statistics in the nearby LF, (iii) 
the use of inaccurate mass- magnitude relationships in the M-dwarf regime. 
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Integrations of this IMF below 1 M Q and of the present-day MF for main sequence 
stars above 1 M Q , plus the contribution of stellar remnants and red giants, yield the 
presently most reliable determination of the disk stellar density and its contribution to the 
Galactic disk mass budget. 

The author is indebted to Pavel Kroupa and Neill Reid for providing their luminosity 
functions and for various discussions. The author is most grateful to the Astronomy 
Department and to the Miller Institute of the University of Berkeley, where most of this 
work was conducted, for a visiting professor position and for their warm hospitality. 
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Table 1: Parameters of the mass functions 



£(m) (eqn.(6),IMFl) f(log m) (eqn.(7),IMF2) f(m) (eqn.(8),IMF3) 



0.141 pc" 3 3.0 M" 1 pc~ 3 

3.3 



A 


0.019 M" 1 pc~ 3 


a 




m < 1.0 M 


1.55 


m > 1.0 M 


2.70 


m {M & ) 




P 




a 





0.1 716.4 

0.25 
0.627 



(m) ( M ) 0.44 0.45 0.55 
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Table 2: Disk present-day stellar density. 





scale height 


n* 


fh 


S 




(pc) 


(pc~ 3 ) 


(M pc" 3 ) 


(M pc- 2 ) 


Thin disk: 










LMS (< l.OM ) 


250 


0.12± 0.02 


3.10 ±0.3 x 10~ 2 


15.5 ±2 


MS stars > 1.0 M 


Scalo'86 


0.43 x 10~ 2 


0.6 x 10~ 2 


2.3 


WD + NS + RG 


250 


0.7 x 10~ 2 


0.6 x 10~ 2 


2.8 


all stars 




0.13 ±0.02 


4.30 ±0.3 x 10~ 2 


20.6 ±2 


Thick disk: 


760-1000 








all stars 






^0.22 ±0.02 x 10~ 2 


^3.2-4.3 


Total 




0.13 ±0.02 


4.50 ±0.3 x 10~ 2 


24.4 ±2.5 
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FIGURE CAPTIONS 



Fig. 1. — Mass functions log[£(logm)] in 10 3 pc 3 derived from the LFs Qv 52p c (solid line), 
$Vs P c (dashed line) and §k 6pc (dotted line), with the BCAH (1998) m-M v and m-M K 
relationships. 



Fig. 2. — Same as Fig. [I] but calculated with the Delfosse et al. (2000) m-My relationship 
for m < 0.6 M & and from the corrected <&k 8 c ( see text). 



Fig. 3. — Same as Fig. ||for the MFs derived from $y 52pc (solid line) and §v apc (dashed line). 
The curves display the 3 different functional forms: power-law (IMF1; solid), log-normal 
(IMF2; dash-dot) and exponential (IMF3; dash) (see text and Table 1), all normalized at 
0.8 M Q . The Kroupa (2001) IMF is shown by the dotted line, as a series of power-laws. 



Fig. 4. — Same as Fig. [3] over the mass-range 0.01 M & < m < 60 M Q . 



Fig. 5. — Same as Fig. ^| in a linear scale £(m) over the mass-range 0.1 M < m < 0.6 M . 
The dotted line corresponds to a Salpeter IMF with the same normalization. 
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